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Explaining cosmic inflation by the effective dynamics of an SU(Nc) pure 
Q_i| gauge theory of scale A ~ 10 -6 M p is lacking an explanation of the (Gaussian) 

O • spatial curvature perturbations needed to seed the formation of large-scale 

structure after inflation. In this work it is demonstrated how fundamentally 
charged fermions of mass mp ~ 10 -12 A, whose approximate chiral symmetry 
' is spontaneously broken during inflation, can cure this shortcoming. The as- 

^ ■ sociated gas of weakly interacting pseudo Nambu-Goldstone bosons (PNGB) 

undergoes Bose-Einstein (BE) condensation well before the end of inflation. 
This causes the occurence of condensed light scalar fields effectively acting 
as a curvaton. Fermions may also be charged under a gauge group G with 
a weak coupling g. Since fermions charged under SU(Nc) are confined af- 
ter inflation the decay of Nambu-Goldstone bosons, which reside in the con- 
densates and in the PNGB radiation generated at reheating, is mainly into 
fermions solely charged under G. The associated decay rate V is estimated us- 
ing PC AC and large Nc counting. PNGB decay takes place during radiation 
domination after cosmological scales have entered the horizon. Neglecting the 
effects of the spontaneous breaking of G induced by the BE condensation, it is 
demonstrated that the observed spectrum of spatial curvature perturbations 
P^ 2 ~ 5 x 1CP 5 is compatible with g ~ 1CP 3 , a ratio r^ ec ~ 10 -4 of curvaton 
to radiation energy at PNGB decay, and a curvaton effective equation of state 

Pcurv — 0.9/9 curtI . 



1 Introduction 



Cosmic inflation is an attractive approach to the horizon, flatness, monopole, and 
homogeneity problem pQ. It is of primary interest to find an answer to the question 
how large-scale structure is seeded and how the anisotropies of the cosmic microwave 
background at large angular resolution come about, see j2] and references therein. 
More recently, string inspired alternatives to inflation have been proposed where 
cosmology is driven by the dilaton field [3] or by a radion field [3] which determines 
the distance between two parallel and four dimensional boundary branes of a five di- 
mensional bulk. In their minimal versions both models, however, have been claimed 
to not be able to generate a scale-invariant spectrum of adiabatic fluctuations E] ■ 

The explanation of inflation as well as the description of adiabatically gener- 
ated primordial curvature perturbations by the real-time dynamics of one or more 
slowly rolling, real, and gravitationally minimally coupled scalar fields is an effective 
approach. Inflationary models of this type are hardly constrained by observation, 
and hence a large amount of arbitraryness is inherent to the model building. To 
attenuate this unsatisfactory fact and to build on a solid principle it was proposed in 
Ref. that the strong dynamics driving very early comology is due to an SU (Nc) 
gauge symmetry with large Nc- 

Effectively, the thermal behavior of the pure gauge theory is modelled by quasi- 
particle excitations over a nontrivial classical background. An adjoint Higgs model 
at very high temperatures is smoothly connected to a dynamical Abelian projection 
at high temperatures which, in turn, reduces to a Z^ c symmetric model for Nc — 1 
complex scalar fields at T < A. jSj. The latter model has a vanishing vacuum 
energy and therefore terminates inflation. Here A denotes the dynamical scale of 
the SU(Nc) theory which also sets the scale for inflation. In this paper we focus 
on the Abelian and the center phase. A finite modulus \<pk\ < A embodies the 
T dependent condensation of magnetic monopoles whereas \<f>k\ ~ A describes the 
condensation of center vortices. Qualitatively, this picture is supported by lattice 
simulations of an SU(2) gauge theory at finite temperature jHj. The potential V 
of our effective description is constructed to admit topologically nontrivial, periodic 
Bogomoln'yi-Prasad-Sommerfield (BPS) 10 j saturated solutions to the Euclidean 
0- field equations along a compact time dimension of length /3 = 1/T [7| lllj. 

The effective gauge theory description implies radiation domination for Mp > 
T > A which runs into vacuum domination for T ~ A The latter regime is 
terminated by a transition to the phase of condensed center vortices. In the absence 
of fermions no (pseudo) Nambu-Goldstone moduli are produced in this transition 
since Z^ c is a discrete symmetry. Heavy vector bosons, which survive inflation, 
are stable and possibly contribute to today's cold dark matter. In the presence of 
fundamentally charged, low-mass fermions a kinematically possible, postinflationary 
decay of vector bosons into Zjv c -charged chiral pseudo Nambu-Goldstone bosons 
(PNGB) should be strongly suppressed by higher dimensional operators. 

Since the monopole condensate fluctuates with a mass ~ A ^> Hi during 
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inflation, denoting the (quasi-constant) Hubble parameter, the usual mechanism 
of an adiabatic generation of spatial curvature perturbations fails in the case of 
pure gauge dynamics [7j. To propose an alternative mechanism for the generation 
of spatial curvature perturbations in the framework of an effective gauge theory 
description of inflation is the purpose of the present paper. The scenario relies on 
the curvaton mechanism which was discovered in Ref. ^3] and re-addressed recently 
in Ref. j!4j . In its minimal version this mechanism operates with a gravitationally 
minimally coupled, light scalar field - the curvaton - which is decoupled from the 
dynamics driving inflation. During inflation the quantum fluctuation of the curvaton 
convert into isocurvature perturbations on super-horizon scales. Upon horizon-entry 
of cosmological scales during the radiation dominated epoch after inflation isocur- 
vature perturbations get converted into spatial curvature perturbations due to an 
increasing ratio of curvaton to radiation energy density. Following the epoch of cur- 
vaton decay into radiation the spatial curvature perturbation remains constant in 
time and is responsible for the formation of large-scale structure [2]. In our gauge 
theory framework it is suggestive that the existence of a curvaton has to do with 
the dynamical breakdown of a continuous, global symmetry due to strong gauge 
interactions close to the transition to the confining phase. This symmetry is present 
if we allow for fundamentally charged, chiral fermions to enter the stage. In this 
case, the spontaneous breakdown of chiral symmetry during inflation releases a gas 
of chiral PNGB which undergoes Bose-Einstein (BE) condensation after a sufficient 
number of e-foldings. The resulting condensate acts as a curvaton and decays into 
fermions solely charged under a weakly coupled gauge group G. 

The outline of the paper is as follows: In Section 2 we introduce our model for 
the relevant fields during inflation, discuss BE condensation of chiral PNGB in an 
expanding universe, and use a simplified, effective description of the condensate to 
estimate the mass of fundamental fermions from the requirement of scale invariance 
of the Gaussian isocurvature perturbations generated during inflation. The postin- 
flationary generation of curvature perturbations and the subsequent decay of PNGB 
are discussed in Section 3, and the consistency of various assumptions is checked. 
There is a short summary in Section 4. 

2 Chiral symmetry breaking and BE condensa- 
tion during inflation 

2.1 Set-up and general considerations 

In this paper we consider the dynamics of matter that is not directly related to the 
Standard Model of particle physics and its supersymmetric extensions. 

We consider Np flavors of asymptotically free fermions of degenerate mass mp 
which transform fundamentally under an SU(Nc) gauge theory of scale A. At one 
loop accuracy we have Nc > tt Np. Assuming nip <C A, there is an approximate 
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chiral symmetry U(Np). For this symmetry group to possess nontrivial, nonanoma- 
lous subgroups we impose N F > 1. To have a sensible l/N c counting we assume 
Np <C Nq- The fermions are taken to be fundamentally charged under an additional, 
much smaller gauge group G inducing weakly coupled dynamics with a gauge cou- 
pling constant g < l 1 . The total gauge group therefore is SU(Nc) x G. We will 
refer to the fermions characterized above as F-fermions. Assuming that the vacuum 
energy of the theory is dominated by that of the pure glue sector, one obtains an 
effective groundstate description as in [7|E]- The magnetic coupling constants of 
the SU(N C ) gauge bosons to the monopole condensate, which is described by a 
single, complex mean field 0, are determined by the root structure of SU(N C ) [T2"] . 
For simplicity we assume all magnetic couplings and hence all gauge boson masses 
to be equal. 

In addition, K flavors of fermions with a mass comparable to are assumed. 
These fermions be fundamentally charged solely under G, and in the following we will 
refer to them as .ff -fermions. Einstein gravity is assumed to be minimally coupled 
to this matter. 

For large Nq the critical temperature for the transition to the confining Z^ c 
symmetric phase is T ~ A- [7j . If we assume that charge parity C is a good symmetry 
in the deconfining phase then close to the transition the ratio R of radiation pressure 
to vacuum pressure is given as 



Eq. (0) derives from the facts that the gauge bosons of the dynamical Abelian pro- 
jection decouple close to the transition to the center phase jH| and that there are 
A(NpNc + Xdim[fund(Gr)]) chiral fermions and 2 dim[adj(G)] massless gauge bosons 
of the gauge group G. 

Hence the ground state dominates the cosmic evolution (quasi-exponential ex- 
pansion, quasi de Sitter geometry) close to the phase transition if K and Np are not 
too large. For example, N c = 34, G = SU(2), and K = Np = 3 yields R ~ 3.6 %. 

Quenched QCD lattice simulations seem to indicate that the critical tempera- 
tures for chiral restoration is slightly higher than the deconfinement transition [T7] . 
Including dynamical fermions and appealing to the method of imaginary chemical 
potentials it was found recently in Ref. ^Hj that Nc = 2,3 QCD has no confinement 
phase transition at all but rather a level crossing. This is consistent with ^1] where 
the nonexistence of a tachyonic regime in the model effective potential was observed 
for Nc < 8. However, the classical arguments of Ref. [TT] are only applicable in the 
limit of large Nc, and therefore this consistency at small Nc may be pure coinci- 
dence. We nevertheless take the results of J7j as a motivation to assume that in the 
framework of Ref. [7j, where Nc = 34 was found from the condition that inflation 
is terminated after 60 e-foldings, the chiral transition occurs during the de Sitter 
stage. This transition naturally generates a radiation of PNGB. 

1 The consistency of this assumption will be shown in Section 4.2. 
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Our intuition about inflationary cosmology seems to forbid the condensation of 
PNGB since the exponential expansion of the universe tears particles apart and 
allows for causal contact only within a quasi-constant horizon. This, however, may 
not be a valid concern: The total particle number as measured at the instant of 
emission is only conserved on super-horizon scales. A particle being emitted at some 
instant to is m a momentum eigenstate normalized to the Hubble volume. At some 
later time t\ > t this volume has expanded by a factor exp[3H(ti — t Q )] decreasing 
the contribution of the particle to the total density by a factor exp[— 3H (ti — to)]. 
Since the position of a particle is irrelevant for the formation of condensates we may 
still think in terms of conserved particle numbers. Therefore, the relevant patch 
of the inflating universe, which coincides with the observable universe at emission, 
may be associated with a box of exponentially increasing volume containing a gas of 
exponentially decreasing temperature. Bose-Einstein (BE) condensation is a pure 
quantum effect which needs no interaction betweeen the participating particles JH] ■ 
Hence a comparison of expansion and interaction rate again is irrelevant for the 
formation of condensates. Therefore, BE condensation may be expected to take 
place after a sufficiently long time during the de Sitter epoch. 

The time scale for the duration of the chiral transition is set by / PNGB ~ A. Since 
A ~ ^jr-Hi ^> Hi this transition takes place on a much smaller time scale than the 
inflationary expansion of the universe, and hence the chiral transition may be viewed 
as instantaneous. We assume that the spectrum of PNGB generated at the instant 
of emission is of the Bose-Einstein type with temperature T ~ A and chemical 
potentials /z" . In a given gauge the correspond to charges within the horizon 
arising from the inhomogeneities in the charge densities n a which are defined as the 
zero components of G gauge covariantly conserved currents j". Because the coupling 
g is small and since the chiral interaction strength between PNGB is parametrically 
counted in powers of 47i -^ NGB (/ PNG b ~ A) and at the instant of PNGB emission 
T ~ A- the treatment of the PNGB as a Bose gas of noninteracting particles is 
justified. 



2.2 BE condensation in an expanding universe 

We have to distinguish two cases. A chiral PNGB particle is either (a) charged or (b) 
neutral under the group G. Condensation of charged PNGB can only occur if there 
are, in a given gauge, nonvanishing charge densities n a within the de Sitter horizon 
being produced during the chiral transition. Here a runs through the representation 
of G which acts on the charged PNGB being made of a pair of an F-fermion and 
F-antifermion. There are Np species of PNGB for each G charge. We do not carry 
this factor along in the following discussion. Moreover, we will omit the label a in 
the following since for generic gauge groups many of the n a can be gauge-rotated 
away. 

Let us first consider ). A G charge Q e within the inflationary horizon at 

the instant of the chiral transition is associated with a chemical potential fi e . We 
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have 



— dpp 2 x 



ex P[f(vP 2 + m PNGB - Me)] - 1 exp[f ( Jp 2 + m2 NGB + fi e )] - 1 



•(2) 



Let x denote the number of e-foldings the scale factor a undergoes from the instant 
of the chiral transition to the first occurence of the charged BEC. Then the charge 
density n e relates to the charge density n at the onset of condensation as 

n Q = exp[— 3x] n e . (3) 

In an ideal, nonrelativistic, and bosonic quantum gas of identical particles the charge 
density no, in turn, determines the critical temperature To for the onset of BE 
condensation as [IH] 

n = 2.612 x(^^) 3/2 . (4) 
Combining Eqs. (0),©, and (JJJ), we derive 



T ° = £b eXph2x] (2.612x2* Jo" ^ 



1 



ex P[x(vP 2 + m PNGB - »e)] ~ 1 exp[^ Up 2 + m2 NGB + fi e )] - 1 



2/3 



(5) 



Eq. (jSJ) relates the critical temperature T to the chemical potential at emission /i e 
and the number of e-foldings x needed for the onset of condensation. 

After condensation there must be sufficient inflation to generate classical, Gaus- 
sian isocurvature fluctuations of the BEC's which are converted to curvature per- 
turbations after inflation (see Section 3.1). Therefore, BE condensation must be 
completed within a number of e-foldings considerably smaller than the total of 60 
or so inflationary e-foldings required to solve the horizon problem. 

It is the finiteness of the scale m PNGB that makes it likely to cause a finite charge 
density n e during the chiral transition. Thus we assume /x e to be comparable to 
m PNGB . Moreover, we assume that m PNGB also sets the scale for T . There is no strict 
reason why this should be so. However, m PNGB = /x e = T leads to an acceptable 
number of e-foldings for the completion of BE condensation if we use A ~ 10 13 GeV 
as in Ref. [7j. 

Let us make this more explicit. As a result of Section 2.4 we take m PNGB = 10~ 6 A. 
Inserting m PNGB = \i e = T , A = 10~ 6 M p , m PNGB = 1CT 6 A into Eq. (JHJ) and varying 
/j, e between \i e = 0.1m PNGB and \i e = 0.9m PNGB , the number of e-foldings x needed 
for the onset of BE condensation varies between x = 7.5 and x = 8.2. Hence, 
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there is a good insensitivity of x with respect to variations in /i e . Note that with 
~ TH ~ ^ 12 M p we satisfy the condition Hi < 10 5 M p which is inferred from 
the cosmic microwave background limit on gravitational waves. Here M p denotes the 
Planck mass M p ~ 10 19 GeV in contrast to a reduced Planck mass Mp introduced 
below. 

To estimate the number of e-foldings x c needed to complete condensation we 
consider the following ratio as it was derived in Ref. ^J] 



ntot \T J 

Here, n tot , n c denote the total density of PNGB and the density fraction of particles 
residing in the condensate, respectively. We decompose n tot (x) = n c (x) + n(x), 
where n(x) refers to the fractional density of particles in the gas phase after x e- 
foldings counted from the onset of condensation. For T < m PNGB temperature scales 
nonrelativistically, that is T = T exp[— 2x]. Substituting this into Eq. (JHJ), we obtain 

n(x) = n c (x) (exp[3x] — l) -1 . (7) 

Since particles in the condensate have no spatial momentum their density n c does 
not scale with x as a direct consequence of the expansion of the universe. Because 
of Eq. (0), which describes the T dependence of the condensate fraction, there is, 
however, an indirect scaling. This can be seen as follows. 

Changing x — > x + dx, the variation dn is due to gas particles going into the 
condensate on the one hand and due to a decrease of gas density because of de Sitter 
expansion on the other hand. Hence, we have 

dn = —(n' c + 3n)dx , (8) 

where n' c = ™. Replacing x — > x + dx in Eq. (f7j) . using Eq. and comparing 
coefficients at dx, we obtain 

expl-6,] 
1 — exp[— 3xJ 

The solution to Eq. (jH)) reads 

n c (x) = n C)OQ exp [exp[— 3x}} (1 — exp[— 3x]) . (10) 

Eqs. ©)© can also be solved for n c ,dn c , respectively. Repeating the above steps 
in this case yields 

n' = -3 eX P[- 3 "U , (11) 

l + exp[3x] v ; 

The solution to Eq. (fTTj) . subject to the initial condition n(0) = n , is 

n(x) = — exp [exp[— 3x] — 3x] . (12) 

e 
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Figure 1: Density of particles, which reside in the condensate, as a function of the 
number of e-foldings x as counted from the onset of condensation. 



From Eq. (JJJ) we finally obtain 



n r 



n 



(13) 



Fig. 1 shows a plot of n c (x)/n CtOQ . Saturation of the condensate density n c is reached 
within a small fraction of an e-folding at a value comparable to n^. Hence, given an 
initial chemical potential the number of e-foldings required for a complete condensa- 
tion of PNGB particles is practically determined by the critical temperature To for 
the onset of the condensation. This shows that the assumptions m PNGB ~ T ~ \i e 
are consistent with the requirement that BE condensation be completed after a 
number of e-foldings which is considerably smaller than 60. From Eqs. P)). (jl3j) it is 
seen that the asymptotic condensate density is 



n,. 



171, ~ 
"'PNGB 



0.11/ G m 



3 

PNGB ) 



(14) 



where the factor /g is determined by the gauge group G and the result of PNGB 
emission. For example, if G — SU(2) then charged PNGB, composed of pairs 
of F and F fermions, would live in the adjoint representation. Assuming equal 
chemical potentials ^ = m PNGB , (a = 1, 2, 3) for each charge in the gauge we start 
with would yield fsu(2) — v3. The condensates of charged PNGB would act as an 
adjoint Higgs field 4> a , (a = 1, 2, 3) and break SU (2) spontaneously. We will address 
this concern in Section 3.3. To obtain the total number density of PNGB residing 
in the condensates the expression in Eq. ()14|) has to be multiplied by Np. 
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Let us now consider case (b). Here we have 



T 



2n 



exp[— 2x] 



2.612 x 2tt 2 



dpp 2 



exp[-fVF 



mi 



2/3 



We assume again that BE condensation sets in at T 



~ m 



PNGB • 



(15) 

The corresponding 



number of e-foldings is ~ 12. As for the rapidity of the process of complete conden- 
sation we expect a similar behavior as for case (a). Setting fa = 1, the asymptotic 
condensate density is again N F times the value estimated in Eq. (JHJ). 



2.3 Simplest field theoretic description of condensates 

Similarily to Ref. [201, where the simplest curvaton model is discussed, we describe 
the flavor-sum over BEC's of PNGB in terms of one effective, condensed real scalar 
field a* which is one of the minima of some potential V{a). Since, according to what 
we assumed in the last section, there is no scale hierarchy in the condensates, the 
contribution to the potential energy of a* is ~ (c*) 4 - Subdominant fluctuations 5a 
give a contribution |m 2 5a 2 to lowest-order, where m 2 a = d%.V\ a=a *. Hence we have 

V^{a*Y + \ml5a\ (16) 

To estimate the value of a* we have to know that the transition to the condensed 
phase has no latent heat, that is, there is no discontinuity in the energy density [T9"] . 
Therefore, the energy density p for a single species of PNGB is given as 2 

P = n c,OQ m PNGB • (17) 

Combining Eq. (fT7J) and Eq. (fTl]) yields 

p = 0.11/ G m 4 NGB . (18) 

Taking into acount flavor multiplicity and the group structure of G, the factor 0.11 
in Eq. (J18)) converts into a number of order unity such that the total energy density 
Pbec is given as 

Pbec ~ mp NG B . (19) 
Combining this with Eq. ()16j) . we have 

a* ~ m PNGB . (20) 

This can be converted into 

a* ~ yjAm F . (21) 

2 Because of the existence of a condensate of charged PNGB and consequently the existence of 
a finite G held strength there is a small contribution to the energy density from the gauge held 
energy which we omit here. 
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when invoking large Nc counting and the Gell-Mann-Oakes-Renner (GOR) relation 

m PNGB /pngb — m F (qq) , (22) 

where (qq) denotes the order parameter for spontaneous chiral symmetry breaking. 
More explicitly, Eq. (|2"Tj) can be obtained as follows. On the one hand large Nq 
counting implies that f PNGB ~ V^c A. Here A oc N c indicates the asymptotically 
constant in Nc scale parameter of the SU {Nc) gauge dynamics with coupling con- 
stant g s /y/Nc |21j. On the other hand we have — (qq) ~ A^A 3 . Since mp does 
not scale with Nc we hence conclude that m PNGB ~ y/mp~A which links Eq. (|2U|) and 
Eq. (ED). 



2.4 Estimate of fermion mass from scale invariance of Gaus- 
sian isocurvature perturbations 

We now discuss the generation of Gaussian isocurvature perturbations due to quan- 
tum fluctuations of the BEC of PNGB during inflation with the simplifying assump- 
tion of an effective description in terms of a single, real curvaton field a. The case 
of non-Gaussianity is not investigated in this paper. 

For the generation of Gaussian isocurvature fluctuations during inflation the 
typical perturbation 5a = Hi/2ir acrued in a Hubble time must be small as compared 
to the mean-field a* as it is given in Eq. (|2"T|) . Therefore, we demand 

4?rV 2 2 Ml m F 

^-4.— (23) 

where Mp denotes a reduced Planck mass defined as H 2 = Mp 2 p. It is natural to 
assume that the mass m a of the fluctuations 5a is comparable to m PNGB . Hence, 
for 5a to be the dominant cause of density perturbations and for the generation of 
a scale invariant spectrum of classical perturbations we have mp NGB ~ m 2 H 2 . 
This is equivalent to 

A 2 A 

-— 77- > 1 . 24 

Mp mp 

Hence, Gaussianity and scale invariance of the spectrum imply 

A 2 m F 1 A 2 , , 

»-^»7-7T7 2 -- (25) 



Mp A 4vr 2 M P 
In particular, assuming Mp/A 2 ~ 10 11 as in Ref. [7], we have 



KT 13 < ^ < 10 -11 , (26) 

which is a tight constraint on mp. Taking ^ ~ 10 -12 , corresponding to m F ~ 
10 GeV, yields the value m PNGB ~ 10~ 6 A which is used in Section 2.2. 
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During inflation the Hubble paramter is given as Hi ~ jg-A. From Eq. flUJ) it 
then follows that 



This is an extremely small number. Therefore, fluctuations of a during inflation 
solely induce isocurvature perturbations. 

3 Curvature perturbations and the decay of PNGB 
3.1 The curvaton mechanism 

In this section we briefly review the curvaton mechanism as it was proposed in 
Ref. [T3J and reconsidered by Lyth and Wands in Ref. [H| (see also [22] )■ We also 
propose a modification concerning the nature of the curvaton matter after inflation. 

Most economically the generation of primordial curvature perturbations is at- 
tributed to the quantum fluctuations of the inflaton field during inflation. In the 
standard approach fluctuations are adiabatically converted into classical, Gaussian 
curvature perturbations several Hubble times after horizon-exit. Upon horizon-entry 
after inflation these curvature perturbations seed the formation of large-scale struc- 
ture. 

Believing in this simple picture, one ignores attractive models of inflation where 
the curvature perturbations cannot be generated in an adiabatic way. A generic 
alternative was discussed in Ref. and applications in more specific models were 
carried out in Refs. [2B] . The curvaton mechanism is based on the idea that primor- 
dial curvature perturbations do not have to be created directly during inflation but 
can be switched on by the cosmic evolution in a postinflationary, radiation dom- 
inated universe. Since in the framework of gauged inflation we have in mind an 
inflation scale A ~ 10 13 GeV [7] the effects of gravitational waves on the generation 
of curvature perturbations are neglected in this paper. 

Let us now turn to an application of the curvaton mechanism to a fluctuating 
BEC of chiral PNGB. A practically scale invariant spectrum of Gaussian isocur- 
vature perturbations, generated by a cosmologically irrelevant light field during 
inflation, may induce primordial curvature perturbations in the course of the cos- 
mic evolution following inflation. How isocurvature perturbations are generated on 
super-horizon scales during inflation was already discussed in Section 3.4. Here we 
address the postinflationary conversion into curvature perturbations. According to 
(|TB|) small fluctuations 5a <C a* satisfy 



Coherent oscillations in 5a start if H becomes of the order m PNGB . In a radiation 
dominated universe this happens at 




(27) 



5a + 3H5a + mi NrK 5a = . 




a 



mass 



(29) 



^end 

\fm mGB M P ' 
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where a en d denotes the scale factor at the end of inflation. With Mp/A ~ 0.35 x 10' 
and m PNGB ~ 10~ 6 A Eq. (J2HD yields 



"mass 



1.7. (30) 

(lend 

Hence, 5a starts to oscillate right after inflation. Since in Eq. ()16j) the mean field part 
and the coherent oscillations 5a act like a cosmological constant and nonrelativistic 
matter, respectively, we assume a mean equation of state for curvatonic matter 
which interpolates between pure matter (p = 0) and vacuum (p = —p) 

Pcurv = ^Pcurv j (0 ^ AC ^ 1) . (31) 

Admittedly, this corresponds to a rather rough treatment of the curvatonic matter 
since the ratio between the nonrelativistic matter part and the condensate part 
does evolve in time. The equation of state Eq. (|3~Tj) thus is merely assumed to be a 
time average over the exact equation of state taken from the end of inflation until 
curvaton decay. 

Continuity implies that 

/ a \ 3(k— l) 

Pcurv Pcurv,end I I ■ \^ ) 

V de-rid / 

Since "baryons" and heavy "mesons" are of mass Nq A, ~ A ^> m PNGB ^T] the 
transition to the confining, Z^ c symmetric phase, which terminates inflation, should 
mainly generate PNGB radiation due to the much larger phase space. If we assume 
the distribution of this radiation to be of the Bose-Einstein type then an extra 
contribution oc <5(p) has to be added to this distribution to take into account the 
BEC of PNGB (curvaton field) already generated during inflation. Consequently, 
the spontaneous breaking of chiral symmetry during inflation contributes to the 
initial conditions at reheating which were postulated in Ref. |24j . 

Immediately after inflation the universe is dominated by PNGB radiation. We 
will show in Section 4.3 that this situation is maintained until the decay of the curva- 
ton. Radiation as well as curvatonic matter separately satisfy continuity equations 

P/3 = —3H(pp + pp) ((3 = rad, curv) . (33) 

Each fluid generates curvature perturbations £^ which are defined on slices of uniform 
P/3 and which are time independent on super-horizon scales. From Eq. (}3*3*j) we have 

6, = -*^ (34) 

P/3 3 pp + Pf3 

The total curvature perturbation £ is time dependent and reads [22] 

/-/.\ Prad(t)£rad + Pcurv (t)£curv 1 r(t) 5p 

curv r o _ \ 

Prad{t) +Prad{t) + Pcurv 

(*) + 

Pcurv 

(t) 3 4 + r(t)(l Ac) Pcurv 
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where r(t) = Pc " r "ff , and the density contrast has been neglected. If we assume 



§£c_ 



a value of k not far from unity then, according to Eq. JED) , the density contrast 

; pcurv 

is given as 

$ Pcurv 1 / -^i 



(36) 

Pcurv 2 \ 27r77l PNGB / 

The consistency of this assumption will be shown in Section 4.3. The curvaton 
decays at a rate V due to a decay of PNGB into A"-fermions (see Section 4.2). 
After complete decay, which we assume happens instantaneously at H = T, £ is 
constant in time. Therefore, the prediction for the primordial spectrum of curvature 
perturbations (the Bardeen parameter) is 



,1/2 _ 1 r dec / H, 

6 4 + r dec (l — k) \27rm 



p e ~ £ t-^t^ — - ( o_l ) • ( 37 ) 



In a radiation dominated universe r dec calculates as 

.Odec/ V A / \a dec 



r dec ~ n I ~ — — . (38) 

\a d ecJ V A / \a dec J 

The measured value of the Bardeen parameter is P^ 2 ~ 5 x 1CT 5 . 



3.2 Decay channels 

The "hadronic" modes not being PNGB, which exist after inflation, are expected 
to possess mass of order A, AcA ^> m PNGB [IT] . Therefore, the decay of PNGB 
into these modes is kinematically impossible even though these "hadrons" may be 
charged under the weakly coupled G sector. The only possible decays of PNGB are: 
(a) decay of charged and neutral PNGB into a pair of i^-fermions and (b) anomaly 
mediated decay of neutral PNGB into a pair of G gauge bosons. 

Let us first consider case (a). We do not consider a spontaneous breakdown of 
G gauge symmetry by the presence of charged BEC of PNGB. Using the results 
of the next section it will be argued below that this assumption can be rendered 
consistent. Therefore, the corresponding gauge fields are taken to be massless. Since 
we do not want to make any specific assumptions about the structure of G we treat 
group factors to be of order unity - in accord with the general level of precision 
intended in this paper. Since K is assumed to be much smaller than A^ we neglect 
flavor degeneracies. 

We assume F-fermion-gauge-field vertices to be parity conserving in the under- 
lying SU(Nc) theory. Therefore, the PNGB only couples to the axial current jjj. 
As usual, this coupling is parametrized by a decay constant / PNG b as 

(0|4|PNGB,p) = z/ PNGBP/i . (39) 

Diagrammatically, the decay of a PNGB into a pair of A'-fermions is shown in Fig. 2. 
For the decay of PNGB into A'-fermions we distinguish two extreme cases. 
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Figure 2: Amplitudes for a charged PNGB decaying into (i) a pair of i^-fermions 
and (ii) into two G-gauge bosons. 



1) There is an exact flavor symmetry, that is, all i^-fermions have equal mass trip, 
but the G-gauge-boson-fermion vertex is of the A-type, gW^j^, and hence it violates 
parity maximally in the Lagrangian. 

2) The gauge-boson-fermion interaction is of the parity conserving V-type, gW^j^, 
but there is a nontrivial CKM matrix which mixes .fT-flavors of different mass with no 
large hierarchy. For a specific decay we denote the typical mass differences between 
i^-fermion and f^-antifermion by Aroj? and the sum of masses by rap. Moreover, 
we assume the elements of the mixing matrix to be of order unity. 

The corresponding decay rate can be calculated on tree-level in analogy to the pion 
decay rate EH] replacing the fermi coupling as -% — > —f- — and considering a finite 

\/2 m PNGB 

neutrino mass. According to the assumptions made under 1) and 2) instead of a 
(V — v4)-type vertex we separately consider A-type and V-type vertices. 

For case 1) the decay rate Ta can be expressed in terms of the gauge coupling 
g, the mass m PNGB , the decay constant / PNGB , and the mass mp of i^-fermions as 

rA = f& f K( 1 _ i j^\ 1,2 ( 1+i j±_\, (40) 

m PNGB \ m PNGB/ \ m PNGB/ 

For case 2) the expression of the decay rate Yy in terms of the gauge coupling g, 
the mass of PNGB m PNGB , the decay constant f PNGB , and Amjr, Trip reads 



9 4 f^r^ Ami / m % \ ( Ami + m% Am 2 p mp \ . , 

rv = f- ™2 — £ i - i £ £ + — . (4i) 



S.-TT 771° \ 711 z / \ 771 ^ 771 

jii " fc PNGB \ "TNGB/ \ '"'PNGB PNGB 
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Since mp , AmF , mF C 1 we will ignore the small corrections in the brackets 
of Eqs. (|4UI41|) in the following. Because there is no large hierarchy in the masses of 



A'-fermions for case 2) we will treat Amp to be of order mp. 

Using the GOR relation of Eq. (|22|) and large Nc scaling, this yields 



IV,a ~ g 4 VmX ~ — ^ g 4 m PNGB = T . (42) 

Eq. (|42)l expresses the fact that the decay width for PNGB into A'-fermions scales 
linearly with Nc and with the fourth power of a small coupling constant j<1. 
Considering case (b), we can use the result of [20] for anomaly mediated decay 



77 ~ 4tt f 2 ' 1 J 



Appealing to Eq. (j22|) and using / PNGB ~ vAc A, we derive 



4 

T^A^I^yVA. (44) 



This leads to 



T 77 2 m F 



(45) 



Hence, decay into AT-fermions mediated by the coupling of the axial current to the 
PNGB is extremely favoured as compared to the anomaly mediated decay into two 
G gauge bosons. 

Let us come back to the issue of spontaneous G symmetry breaking by charged 
BEG This symmetry breaking induces a mass m\y for the G gauge bosons and a 
mass shift Amp for A'-fermions. We expect and Amp to be of order ga* and 
g 2 cr*, respectively. On the one hand, the fmiteness of m^y affects T through a shift 
m PNGB > m PNGB + m w i n ^ ne G gauge-boson propagator. In Section 4.3 we find that 
values of g smaller than 10~ 3 are allowed suggesting that mw ~ 9&* ~ 9 Wpngb < 
10 _3 m PNGB can be neglected in the calculation of Ty,A- On the other hand, we 
have Amp ~ g 2 a* ~ g 2 m PNGB ~ g 2 ^m F A ~ g 2 10 6 mp. Thus, Am F ~ mjr, and 
hence we may also neglect Amp. For a sudden decay at T ~ 10~ 6 A (see next 
section) a thermal mass ~ gT for G gauge bosons can be neglected. However, 
after thermalization the thermal mass of the released quasi-particle excitations of 
A'-fermions is a factor of ~ 10 3 larger than mp. 



3.3 How large are g, a dec /a end , and 

' dec ' 

In this section we combine the conversion mechanism of Section 4.1 and the above 
approximate determination of the PNGB decay rate T to estimate the gauge coupling 
g, adec/aend, and r^c as functions of k. It is assumed that the running of g with 
temperature T can be neglected. The condition for sudden curvaton decay is H = T 
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Figure 3: The gauge coupling g at PNGB decay as a function of k. The curve 
corresponds to P\ 12 = 5 x 10" 5 , A = 10 13 GeV, m PNGB = 10" 6 A, and N c = 34. 



Assuming radiation domination until decay and appealing to Eq. (|42j) . we 
obtain 

adec 1 J*L - = (46) 

(lend 9 2 V N C V M P™>PNGB 

Substitution of Eq. flUJ) into Eq. flUS)) yields 



3k/2+7/2 „ 3/2«+l/2 

i-n.:b / \ A 3(K_1) . (47) 



_ m PB ™ / ott 

''"dec 



^6«+2 VM P A" ( 



We proceed by substituting Eq. (JTTj) into Eq. (|37jl and by solving the result for 
This yields 



A 2 \ 2 1 



pV2 \2irMpm PNOB J 4 



V(a K+ 2) ^ xl/4 



MpN ( 



c 



(48 



1/2 

where is evaluated at curvaton decay. In Figs. 3, 4, and 5 graphical presenta- 
tions are shown of Eq. PH]). Eq. ()46|) . and Eq. (|4T|). respectively. 

We are now in a position to check whether our assumptions actually are consis- 
tent about (i) PNGB radiation dominance at curvaton decay and (ii) the requirement 
g < 1CT 3 for a negligible mass shift of i^-fermions due to the spontaneous G gauge 
symmetry breaking by BEG Let us first look at (i). On the one hand, the PNGB 
gas can be considered radiation only as long as its temperature T does not fall below 
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adet/aend 



500000 





0.8 0.825 0.85 0.875 0.9 0.9250.950.975 



K 

Figure 4: The scale factor at curvaton decay as a function of k. Parameter values 
are taken as in Fig. 3. 



r dec 
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. 000197392L 



0.2 



0.4 



Figure 5: The ratio of the energy densities of curvaton matter and PNGB radiation 
just before curvaton decay as a function of k. Parameter values are taken as in 
Fig. 3. 
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m PNGB . At curvaton decay we thus demand 



"' /C <^— • (49) 



&end m PNGB 



Since mpNGB ~ 10 Fig. 4 indicates that the inequality ()49|) can be satisfied for 
k > 0.9. On the other hand, PNGB radiation dominates the energy density of 
the universe if r dec -C 1. According to Fig. 5 this is well satisfied for the entire 
indicated range of k values. Therefore, PNGB radiation dominance at curvaton 
decay is possible. 

From Fig. 3 we see that (ii) is met for k < 0.9. Hence, (i) and (ii) together fix k 
to be k ~ 0.9. 

Fig. 5 indicates that the ratio of curvaton to radiation energy increased from 
Tend ~ 10~ 24 at the end of inflation to rd ec ~ 10~ 4 at curvaton decay. This is the 
feature that allows for a conversion of the isocurvature perturbations, imprinted in 
the curvaton field at inflation, into curvature perturbations after a postinflationary 
expansion of ~ 10 6 . 

Due to the extremely small coupling g ~ 10~ 3 cosmology would be dominated by 
a plasma of K fermions and G gauge bosons after curvaton decay. The K fermions 
of mass ~10GeV are wimpy and therefore potential contributers to the cold dark 
matter of today's universe. 



4 Summary 

The purpose of this paper was to propose and test a generation and decay mechanism 
for a light field that can play the role of a curvaton in the framework of gauged 
inflation. The idea is that the spontaneous breaking of a chiral symmetry during 
inflation and the subsequent Bose-Einstein condensation of the associated pseudo 
Nambu-Goldstone bosons dynamically generate such an effective light scalar field. 
Postinflationary decay of this field and its excitations into fermions is mediated by 
weakly coupled gauge dynamics. In this scenario the mass scales are fixed by chiral 
symmetry and the requirement of the Gaussianity of the isocurvature perturbations 
that are generated from the quantum fluctuations of the curvaton during inflation. 
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